      REAL FUNCTION DGAUSL*8(FUNC,N)                                    00010000
      IMPLICIT REAL*8(A-H,O-Z)                                          00020000
CCCCC ---------GAUSSIAN QUADRATURE LAGUERRE 4N POINTS(N=1,6)------------00030000
CCCCC          COPIED FROM SYS1.MSL(DGAUSL) IN 1978                     00040000
      DIMENSION A(6,24),X(6,24)                                         00050000
     1,A1(19),A2(19),A3(19),A4(19),A5(19),A6(19),A7(19),A8(11)          00060000
     2,X1(19),X2(19),X3(19),X4(19),X5(19),X6(19),X7(19),X8(11)          00070000
      EQUIVALENCE (A(1),A1(1)),(A(20),A2(1)),(A(39),A3(1)),(A(58),A4(1))00080000
     1,(A(77),A5(1)),(A(96),A6(1)),(A(115),A7(1)),(A(134),A8(1))        00090000
     2,(X(1),X1(1)),(X(20),X2(1)),(X(39),X3(1)),(X(58),X4(1))           00100000
     3,(X(77),X5(1)),(X(96),X6(1)),(X(115),X7(1)),(X(134),X8(1))        00110000
      DATA  X1     /              0.32254768961939230D+00               00120000
     1,  0.17027963230510100D+00, 0.11572211735802070D+00               00130000
     2,  0.87649410478927840D-01, 0.70539889691988760D-01               00140000
     3,  0.59019852181507980D-01, 0.17457611011583470D+01               00150000
     4,  0.90370177679937990D+00, 0.61175748451513070D+00               00160000
     5,  0.46269632891508080D+00, 0.37212681800161150D+00               00170000
     6,  0.31123914619848370D-00, 0.45366202969211280D+01               00180000
     7,  0.22510866298661310D+01, 0.15126102697764190D+01               00190000
     8,  0.11410577748312270D+01, 0.91658210248327360D+00               00200000
     9,  0.76609690554593660D+00, 0.93950709123011330D+01 /             00210000
      DATA  X2     /              0.42667001702876590D+01               00220000
     1,  0.28337513377435070D+01, 0.21292836450983810D+01               00230000
     2,  0.17073065310283440D+01, 0.14255975908036130D+01               00240000
     3,  0.D0                   , 0.70459054023934660D+01               00250000
     4,  0.45992276394183480D+01, 0.34370866338932070D+01               00260000
     5,  0.27491992553094320D+01, 0.22925620586321900D+01               00270000
     6,  0.D0                   , 0.10758516010180990D+02               00280000
     7,  0.68445254531151770D+01, 0.50780186145497680D+01               00290000
     8,  0.40489253138508870D+01, 0.33707742642089980D+01               00300000
     9,  0.D0                   , 0.15740678641278000D+02 /             00310000
      DATA  X3     /              0.96213168424568660D+01               00320000
     1,  0.70703385350482340D+01, 0.56151749708616170D+01               00330000
     2,  0.46650837034671710D+01, 0.D0                                  00340000
     3,  0.22863131736889260D+02, 0.13006054993306350D+02               00350000
     4,  0.94383143363919380D+01, 0.74590174536710630D+01               00360000
     5,  0.61815351187367650D+01, 0.D0                                  00370000
     6,  0.D0                   , 0.17116855187462260D+02               00380000
     7,  0.12214223368866160D+02, 0.95943928695810970D+01               00390000
     8,  0.79275392471721520D+01, 0.D0                                  00400000
     9,  0.D0                   , 0.22151090379397010D+02 /             00410000
      DATA  X4     /              0.15441527368781620D+02               00420000
     1,  0.12038802546964320D+02, 0.99120980150777060D+01               00430000
     2,  0.D0                   , 0.D0                                  00440000
     3,  0.28487967250984000D+02, 0.19180156856753140D+02               00450000
     4,  0.14814293442630740D+02, 0.12146102711729770D+02               00460000
     5,  0.D0                   , 0.D0                                  00470000
     6,  0.37099121044466920D+02, 0.23515905693991910D+02               00480000
     7,  0.17948895520519380D+02, 0.14642732289596670D+02               00490000
     8,  0.D0                   , 0.D0                                  00500000
     9,  0.D0                   , 0.28578729742882140D+02 /             00510000
      DATA  X5     /              0.21478788240285010D+02               00520000
     1,  0.17417992646508980D+02, 0.D0                                  00530000
     2,  0.D0                   , 0.D0                                  00540000
     3,  0.34583398702286630D+02, 0.25451702793186910D+02               00550000
     4,  0.20491460082616420D+02, 0.D0                                  00560000
     5,  0.D0                   , 0.D0                                  00570000
     6,  0.41940452647688330D+02, 0.29932554631700610D+02               00580000
     7,  0.23887329848169730D+02, 0.D0                                  00590000
     8,  0.D0                   , 0.D0                                  00600000
     9,  0.51701160339543320D+02, 0.35013434240479000D+02 /             00610000
      DATA  X6     /              0.27635937174332710D+02               00620000
     1,  0.D0                   , 0.D0                                  00630000
     2,  0.D0                   , 0.D0                                  00640000
     3,  0.40833057056728570D+02, 0.31776041352374720D+02               00650000
     4,  0.D0                   , 0.D0                                  00660000
     5,  0.D0                   , 0.D0                                  00670000
     6,  0.47619994047346500D+02, 0.36358405801651620D+02               00680000
     7,  0.D0                   , 0.D0                                  00690000
     8,  0.D0                   , 0.D0                                  00700000
     9,  0.55810795750063900D+02, 0.41451720484870770D+02 /             00710000
      DATA  X7     /              0.D0                                  00720000
     1,  0.D0                   , 0.D0                                  00730000
     2,  0.D0                   , 0.66524416525615750D+02               00740000
     3,  0.47153106445156320D+02, 0.D0                                  00750000
     4,  0.D0                   , 0.D0                                  00760000
     5,  0.D0                   , 0.D0                                  00770000
     6,  0.53608574544695060D+02, 0.D0                                  00780000
     7,  0.D0                   , 0.D0                                  00790000
     8,  0.D0                   , 0.D0                                  00800000
     9,  0.61058531447218760D+02, 0.D0                    /             00810000
      DATA  X8     /              0.D0                                  00820000
     1,  0.D0                   , 0.D0                                  00830000
     2,  0.D0                   , 0.69962240035105030D+02               00840000
     3,  0.D0                   , 0.D0                                  00850000
     4,  0.D0                   , 0.D0                                  00860000
     5,  0.D0                   , 0.81498279233948880D+02 /             00870000
      DATA  A1     /              0.60315410434163360D+00               00880000
     1,  0.36918858934163750D+00, 0.26473137105544320D+00               00890000
     2,  0.20615171495780100D+00, 0.16874680185111390D+00               00900000
     3,  0.14281197333478190D+00, 0.35741869243779970D+00               00910000
     4,  0.41878678081434300D+00, 0.37775927587313800D+00               00920000
     5,  0.33105785495088420D+00, 0.29125436200606830D+00               00930000
     6,  0.25877410751742390D+00, 0.38887908515005380D-01               00940000
     7,  0.17579498663717180D+00, 0.24408201131987760D+00               00950000
     8,  0.26579577764421410D+00, 0.26668610286700130D+00               00960000
     9,  0.25880670727286980D+00, 0.53929470556132740D-03 /             00970000
      DATA A2      /              0.33343492261215650D-01               00980000
     1,  0.90449222211680930D-01, 0.13629693429637750D+00               00990000
     2,  0.16600245326950680D+00, 0.18332268897777800D+00               01000000
     3,  0.D0                   , 0.27945362352256730D-02               01010000
     4,  0.20102381154634100D-01, 0.47328928694125220D-01               01020000
     5,  0.74826064668792370D-01, 0.98166272629918890D-01               01030000
     6,  0.D0                   , 0.90765087733582130D-04               01040000
     7,  0.26639735418653160D-02, 0.11299900080339450D-01               01050000
     8,  0.24964417309283220D-01, 0.40732478151408650D-01               01060000
     9,  0.D0                   , 0.84857467162725310D-06 /             01070000
      DATA  A3     /              0.20323159266299940D-03               01080000
     1,  0.18490709435263110D-02, 0.62025508445722370D-02               01090000
     2,  0.13226019405120160D-01, 0.D0                                  01100000
     3,  0.10480011748715100D-08, 0.83650558568197990D-05               01110000
     4,  0.20427191530827850D-03, 0.11449623864769080D-02               01120000
     5,  0.33693490584783040D-02, 0.D0                                  01130000
     6,  0.D0                   , 0.16684938765409100D-06               01140000
     7,  0.14844586873981300D-04, 0.15574177302781200D-03               01150000
     8,  0.67216256409354790D-03, 0.D0                                  01160000
     9,  0.D0                   , 0.13423910305150040D-08 /             01170000
      DATA  A4     /              0.68283193308712000D-06               01180000
     1,  0.15401440865224920D-04, 0.10446121465927520D-03               01190000
     2,  0.D0                   , 0.D0                                  01200000
     3,  0.30616016350350210D-11, 0.18810248410796730D-07               01210000
     4,  0.10864863665179820D-05, 0.12544721977993330D-04               01220000
     5,  0.D0                   , 0.D0                                  01230000
     6,  0.81480774674262420D-15, 0.28623502429738820D-09               01240000
     7,  0.53301209095567150D-07, 0.11513158127372800D-05               01250000
     8,  0.D0                   , 0.D0                                  01260000
     9,  0.D0                   , 0.21270790332241030D-11 /             01270000
      DATA  A5     /              0.17579811790505820D-08               01280000
     1,  0.79608129591336300D-07, 0.D0                                  01290000
     2,  0.D0                   , 0.D0                                  01300000
     3,  0.62979670025178680D-14, 0.37255024025123210D-10               01310000
     4,  0.40728589875500000D-08, 0.D0                                  01320000
     5,  0.D0                   , 0.D0                                  01330000
     6,  0.50504737000355130D-17, 0.47675292515781910D-12               01340000
     7,  0.15070082262925850D-09, 0.D0                                  01350000
     8,  0.D0                   , 0.D0                                  01360000
     9,  0.41614623703328550D-21, 0.33728442433624380D-14 /             01370000
      DATA  A6     /              0.39177365150584510D-11               01380000
     1,  0.D0                   , 0.D0                                  01390000
     2,  0.D0                   , 0.D0                                  01400000
     3,  0.11550143395003990D-16, 0.68941810529580860D-13               01410000
     4,  0.D0                   , 0.D0                                  01420000
     5,  0.D0                   , 0.D0                                  01430000
     6,  0.15395221405823440D-19, 0.78198003824594480D-15               01440000
     7,  0.D0                   , 0.D0                                  01450000
     8,  0.D0                   , 0.D0                                  01460000
     9,  0.52864427255691580D-23, 0.53501888130100380D-17 /             01470000
      DATA  A7     /              0.D0                                  01480000
     1,  0.D0                   , 0.D0                                  01490000
     2,  0.D0                   , 0.16564566124990230D-27               01500000
     3,  0.20105174645555030D-19, 0.D0                                  01510000
     4,  0.D0                   , 0.D0                                  01520000
     5,  0.D0                   , 0.D0                                  01530000
     6,  0.36057658645529590D-22, 0.D0                                  01540000
     7,  0.D0                   , 0.D0                                  01550000
     8,  0.D0                   , 0.D0                                  01560000
     9,  0.24518188458784030D-25, 0.D0                    /             01570000
      DATA  A8     /              0.D0                                  01580000
     1,  0.D0                   , 0.D0                                  01590000
     2,  0.D0                   , 0.40883015936806580D-29               01600000
     3,  0.D0                   , 0.D0                                  01610000
     4,  0.D0                   , 0.D0                                  01620000
     5,  0.D0                   , 0.55753457883283570D-34 /             01630000
      K=N                                                               01640000
      IF(N-1)  600,100,100                                              01650000
  100 IF(N-6) 200,200,600                                               01660000
  200 S=0.D0                                                            01670000
      NN=4*N                                                            01680000
      DO 1 J=1,NN                                                       01690000
    1 S = S+A(N,NN-J+1)*FUNC(X(N,NN-J+1))                               01700000
      DGAUSL=S                                                          01710000
      N=K                                                               01720000
      RETURN                                                            01730000
  600 WRITE(6,700) N                                                    01740000
  700 FORMAT(1H0,'(SUBR. DGAUSL) N=',I3,',  N MUST BE LARGER THAN 0 AND 01750000
     1SMALLER THAN 7.')                                                 01760000
      N=6                                                               01770000
      GO TO 200                                                         01780000
      END                                                               01790000
